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A VECTOR ESTIMATION SYSTEM, METHOD AND ASSOCIATED ENCODER 

FIELD OF THE INVENTION 

5 This invention relates to an encoder and a vector 

estimation system and method for processing a sequence of 
vectors to determine a filtered estimate value of each 
vector. The invention is particularly useful for, but not 
necessarily limited to, determining filtered estimate 
10 values to be encoded by a speech encoder and ' transmitted 
over a communication link. 

BACKGROUND ART 

15 A digital speech communication or storage system 

typically uses a speech encoder to produce a parsimonious 
representation of the speech signal. A corresponding 
decoder is used to generate an approximation to the speech 
signal from that representation. The combination of the 

20 encoder and decoder is known in the art as a speech codec. 
As will be apparent to a person skilled in the art, many 
segments of speech signals contain guasiperiodic waveforms. 
Accordingly, consecutive cycles of quasiperiodic waveforms 
can be considered, and processed, by a speech codec as data 

25 vectors that evolve slowly over time. 

An important element of a speech codec is the way it 
exploits correlation between consecutive cycles of 
quasiperiodic waveforms. Frequently, correlation is 

30 exploited by transmitting a single cycle of the waveform, 
or of a filtered version of the waveform, only once every 
20-30 ms, so that a portion of the data is missing in the 
received signal. In a typical decoder the missing data is 
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determined by interpolating between samples of the 
transmitted cycles . 

In general, the use of interpolation by a speech 
5 decoder to generate data between the transmitted cycles 
only produces an adequate approximation to the speech 
signal if the speech signal really is quasiperiodic , or, 
equivalently, if the vectors representing consecutive 
cycles of the waveform evolve sufficiently slowly. 

10 However, many segments of speech contain noisy signal 
components, and this results in comparatively rapid 
evolution of the waveform cycles. In order for waveform 
interpolation in an encoder to be useful for such signals, 
it is necessary to extract a sufficiently quasiperiodic 

15 component from the noisy signal in the encoder. This 
extracted component may be encoded by transmitting only 
selected cycles and decoded by interpolation in the manner 
described above. The remaining noisy component may also be 
encoded using other appropriate techniques and combined 

20 with the quasiperiodic component in the decoder. 

Linear low pass filtering a sequence of vectors 
representing consecutive cycles of speech in the time 
dimension is well known in the speech coding literature. 

25 The difficulty with this approach is that in order to get 
good separation of the slowly and rapidly evolving 
components, the low pass filter frequency response must 
have a sharp roll-off. This requires a long impulse 
response, which necessitates an undesirably large filter 

30 delay. 

A Kalman filter technique for estimating quasiperiodic 
signal components has been described by Gruber and Todtli 
(IEEE Trans Signal Processing, Vol. 42, No. 3, March 1994, 
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pp 552-562) . However, because this Kalman filter technique 
is based on a linear dynamic system model of a frequency 
domain representation of the signal, it is unnecessarily 
complex. It also assumes that the dynamic system model 
5 parameters (i.e. noise energy and the harmonic signal gain) 
are known. However, when considering speech coding, noise 
energy and the harmonic signal gain parameters are not 
known . 

10 A technique for determining the system parameters 

required in a Kalman filter using an Expectation 
Maximisation algorithm has been described in a more general 
setting by Digalakis et al (IEEE Trans Speech and Audio 
Processing, Vol. 1, No. 4, Oct. 1993, pp 431-442). However, 

15 the technique is iterative, and in the absence of good 
initial estimates may converge slowly. It may also produce 
a result that is not globally optimal. No prior art method 
is known for obtaining good initial estimates. Further, 
this method typically requires a significant amount of 

20 data, over which the unknown parameters are constant. In 
the context of speech coding, where the parameters change 
continuously, rapid estimation is essential, and therefore 
this method of applying the Expectation Maximization 
algorithm needs to be improved. 

25 

Stachurski (PhD Thesis, McGill University, Montreal 
Canada, 1997) proposed a technique for estimating 
quasiperiodic signal components of a speech signal. This 
method involves minimizing a weighted combination of 

30 estimated noise energy and a measure of rate of change in 
the quasiperiodic component. This method is highly complex 
and does not allow the rate of evolution of the 
quasiperiodic component to be specified independently. Nor 
does it allow for an independently varying gain for the 

35 quasiperiodic component. 
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In this specification, including the claims, the terms 
comprises, comprising or similar terms are intended to mean 
a non-exclusive inclusion, such that a method or apparatus 
5 that comprises a list of elements does not include those 
elements solely, but may well include other elements not 
listed. 

SUMMARY OF THE INVENTION 

10 

According to one aspect of the invention there is 
provided a vector estimation system for processing a 
sequence of input vectors, said input vectors each 
comprising a plurality of elements, and said system 
15 comprising: 

a digital filter with a filter vector input for 
receiving said sequence of input vectors and a 
predictor gain input for controlling characteristics 
of said filter, said digital filter also having both a 
20 current slowly evolving filter estimate output and a 

previous slowly evolving filter estimate output, said 
current slowly evolving filter estimate output 
providing a current filtered estimate value of a 
slowly evolving component of said sequence of input 
25 vectors and said previous slowly evolving filter 

estimate output providing a previous filtered estimate 
value of said slowly evolving component of said 
sequence of input vectors ; and 

a parameter estimator having an estimator vector 
30 input for receiving said sequence of input vectors and 

a previous slowly evolving filter estimate input 
coupled to said previous slowly evolving filter 
estimate output, said parameter estimator further 
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includes a predictor gain output coupled to said 
predictor gain input, 

wherein when said vector estimation system 
receives a current input vector that is one of said 
5 sequence of said input vectors, said parameter 

estimator provides a current predictor gain value at 
said predictor gain output thereby modifying said 
current filtered estimate value at said current slowly 
evolving filter estimate output, said current 
10 predictor gain value being dependent upon both said 

previous filtered estimate value received at said 
slowly evolving filter estimate input and said current 
input vector received at said estimator vector input. 

15 Suitably, said parameter estimator is characterised by 

said current predictor gain value being dependent upon both 
a sequence of previous said input vectors and a sequence of 
filtered estimate values provided by successive preceding 
values of said previous filtered estimate value. 

20 

Preferably, said current predictor gain value is 
determined by said parameter estimator calculating the 
following : 

25 (Y n T • Xf,^) / (Xf, n -1 T • X f,n-1> 

wherein, y n T is the transpose of said current input vector 
y n that is an nth one of said sequence of input vectors; 
and x f n l T is the transpose of the previous filtered 
30 estimate value x £ n _ 1 resulting from a previous input vector 

Preferably, said filter has a predictor error variance 
output and an observation noise variance input, said 
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predictor error variance output providing a current 
predictor error variance value. 

Suitably, said parameter estimator has an observation 
5 noise variance output coupled to said observation noise 
variance input, and a predictor error variance input 
coupled to said predictor error variance output, said 
predictor error variance output providing a current 
predictor error variance value, 

10 wherein when said vector estimation system receives 

said current input vector, said parameter estimator 
provides a current observation noise variance value at said 
observation noise variance output thereby modifying said 
current filtered estimate value at said current slowly 

15 evolving filter estimate output, said current observation 
noise variance value being dependent upon said previous 
filtered estimate value received at said previous slowly 
evolving filter estimate input, said current input vector 
received at said estimator vector input, said current 

20 predictor gain value and said current predictor error 
variance value . 

Suitably, said current observation noise variance 
value is determined by calculating the following: 

25 

&n -<y a - a n- X f,n-l> /N > ~ S P," 

wherein N is a number of elements of said current input 
vector y n ; S p n is the current predictor error variance 
value associated with said current input vector y ; a is 
30 said current predictor gain value; and x f n _ 1 is said 
previous filtered estimate value. 

Preferably, said parameter estimator has an OnsetFlag 
output coupled to an OnsetFlag input of said digital 
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filter, wherein if a signal at said OnsetFlag input is 
above a threshold value dependent upon harmonic energy in 
said current input vector, said previous filtered estimate 
value is set to a filtered estimate value x f n _ 1 . 

Suitably, if said signal at said OnsetFlag input is 
below the threshold value dependent upon harmonic energy in 
said current input vector, said previous filtered estimate 
value is set to a previous input vector y n _ 1 . 

Preferably, the parameter estimator has an unvoiced 
speech module that determines the current input vector's 
harmonic energy content by assessing the current predictor 
gain value and depending upon the current predictor gain 
value the parameter estimator selectively sets the current 
predictor gain value and the current observation noise 
variance value . 

According to another aspect of the invention there is 
20 provided a vector estimation system for processing a 
sequence of input vectors, said input vectors each 
comprising a plurality of elements, and said system 
comprising : 

a digital filter with a filter vector input for 
25 receiving said sequence of input vectors and an 

observation noise variance input for controlling 
characteristics of said filter, said digital filter 
also having a current slowly evolving filter estimate 
output, a predictor error variance output and a 
30 previous slowly evolving filter estimate output, said 

current slowly evolving filter estimate output 
providing a current filtered estimate value of a 
slowly evolving component of said sequence of input 
vectors, said predictor error variance output 
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providing a current predictor error variance value and 
said previous slowly evolving filter estimate output 
providing a previous filtered estimate value of said 
slowly evolving component of said sequence of input 
5 vectors ; and 

a parameter estimator having an estimator vector 
input for receiving said sequence of input vectors and 
a previous slowly evolving filter estimate input 
coupled to said previous slowly evolving filter 

10 estimate output, said parameter estimator further 

includes a observation noise variance output coupled 
to said observation noise variance input, 

wherein when said vector estimation system 
receives a current input vector that is one of said 

15 sequence of said input vectors, said parameter 

estimator provides a current observation noise 
variance value at said observation noise variance 
output thereby modifying said current filtered 
estimate value at said current slowly evolving filter 

20 estimate output, said current observation noise 

variance value being dependent upon said current input 
vector, said current predictor error variance value, 
and said previous filtered estimate value. 

25 Suitably, said current observation noise variance 

value is determined by calculating the following: 

(Y n T • (Y n - x ffn _ 1 ) /N) - £p, n 

T 

wherein y n is the transpose of said current input vector y n 
30 that is an nth one of said sequence of input vectors; N is 
a number of elements of said current input vector y n ; Zp^ 
is the current predictor error variance value associated 
with said current input vector y n ; and x £ n _ 1 is said 
previous filtered estimate value. 
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Preferably, said parameter estimator has an OnsetFlag 
output coupled to an OnsetFlag input of said digital 
filter, wherein if a signal at said OnsetFlag input is 
5 above a threshold value dependent upon harmonic energy in 
said current input vector, said previous filtered estimate 
value is set to a filtered estimate value x f n _ 1 . 

Suitably, if said signal at said OnsetFlag input is 
10 below the threshold value dependent upon harmonic energy in 
said current input vector, said previous filtered estimate 
value is set to a previous input vector y n _ x . 

Preferably, the parameter estimator has an unvoiced 
15 speech module that determines the current input vector's 
harmonic energy content by assessing the current predictor 
gain value and depending upon the current predictor gain 
value the parameter estimator selectively sets the current 
predictor gain value and the current observation noise 
20 variance value . 

Suitably, said digital filter further includes: 
a slowly evolving predicted estimate output providing a 
current predicted estimate value of said sequence of 
25 vectors. The digital filter may also have a process noise 
variance input. 

Suitably, there may be a smoother module having inputs 
coupled respectively to at least two outputs of said 
30 digital filter. 

Preferably, said smoother module has five inputs 
coupled to respective outputs of said filter. Preferably, 
said smoother module has a smoothed estimate output 
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providing a smoothed estimate value of a previous slowly- 
evolving component. 

Suitably, said smoothed estimate output is coupled to 
5 a smoothed estimate input of said parameter estimator. 

According to another aspect of the invention there is 
provided a method for processing a sequence of input 
vectors each comprising a plurality of elements, said 
10 vectors being applied to a vector estimation system having 
a parameter estimator coupled to a digital filter, said 
method comprising the steps of: 

receiving said sequence of input vectors at 
inputs of said filter and said parameter estimator; 
15 determining a current predictor gain value from a 

current input vector that is one of said sequence of 
said input vectors, said determining being effected by 
said parameter estimator ; and 

applying said current predictor gain value to 
20 said digital filter to thereby modify a current 

filtered estimate value provided at an output of said 
digital filter, said current predictor gain value 
being dependent upon a previous filtered estimate 
value from said filter and said current input vector. 

25 

Preferably, said step of determining may be further 
characterised by providing a current observation noise 
variance value and a current predictor error variance value 
from said current input vector, 

30 

Suitably, said step of applying may be further 
characterised by said filter receiving current observation 
noise variance value thereby modifying said current 
filtered estimate value, said current observation noise 
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variance value being dependent upon said previous filtered 
estimate value, said current input vector, said current 
predictor gain value and said current predictor error 
variance value. 

5 

According to another aspect of the invention there is 
provided a method for processing a sequence of input 
vectors each comprising a plurality of elements, said 
vectors being applied to a vector estimation system having 
10 a parameter estimator coupled to a digital filter, said 
method comprising the steps of: 

receiving said sequence of input vectors at 
inputs of said filter and said parameter estimator; 

determining a current observation noise variance 
15 value from a current input vector that is one of said 

sequence of said input vectors, said determining being 
effected by said parameter estimator; and 

applying said current observation noise variance 
value to said digital filter to thereby modify a 
20 current filtered estimate value provided at an output 

of said digital filter, said current observation noise 
variance value being dependent upon said current input 
vector, a current predictor error variance value from 
said filter, and a previous filtered estimate value 
25 from said filter. 

Preferably, the filter is a Kalman filter. 

According to another aspect of the invention there is 
30 provided an encoder for processing a speech signal, said 
encoder comprising: 

A signal normalization module for processing the 
speech signal to provide a sequence of input vectors each 
comprising a plurality of elements; 
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a digital filter with a filter vector input 
coupled to an output of the signal normalization 
module for receiving said sequence of input vectors, 
the digital filter also having an observation noise 
5 variance input for controlling characteristics of said 

filter, said digital filter also having a current 
slowly evolving filter estimate output, a predictor 
error variance output and a previous slowly evolving 
filter estimate output, said current slowly evolving 

10 filter estimate output providing a current filtered 

estimate value of a slowly evolving component of said 
sequence of input vectors, said predictor error 
variance output providing a current predictor error 
variance value and said previous slowly evolving 

15 filter estimate output providing a previous filtered 

estimate value of said slowly evolving component of 
said sequence of input vectors; 

a parameter estimator having an estimator vector 
input for receiving said sequence of input vectors and 

20 a previous slowly evolving filter estimate input 

coupled to said previous slowly evolving filter 
estimate output, said parameter estimator further 
includes a observation noise variance output coupled 
to said observation noise variance input; and 

25 a slowly evolving component encoder with an input 

coupled to said slowly evolving filter estimate 
output, wherein when said vector estimation system 
receives a current input vector that is one of said 
sequence of said input vectors, said parameter 

30 estimator provides a current observation noise 

variance value at said observation noise variance 
output thereby modifying said current filtered 
estimate value at said current slowly evolving filter 
estimate output, said current observation noise 
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variance value being dependent upon said current input 
vector, said current predictor error variance value, 
and said previous filtered estimate value. 



5 Preferably, the encoder may include an adder module 

with one input coupled said slowly evolving filter estimate 
output and another input coupled to the output of the 
signal normalization module, wherein in use said adder 
subtracts the said current filtered estimate value at the 
10 output of the vector estimation system from at least one of 
the elements of the sequence of input vectors . 



Suitably, an output of the adder module is coupled to 
a rapidly evolving component encoder. 

15 

Suitably, said parameter estimator is characterised by 
said current observation noise variance value being 
dependent upon both a sequence of previous said input 
vectors and a sequence of filtered estimate values provided 
20 by successive preceding values of said previous filtered 
estimate value. 



BRIEF DESCRIPTION OF THE DRAWINGS 



25 In order that the invention may be readily understood 

and put into practical effect, reference will now be made 
to a preferred embodiment as illustrated with reference to 
the accompanying drawings in which: 
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FIG. 1 illustrates a vector estimation system for 
processing a sequence of input vectors in accordance 
with a preferred embodiment of the invention; 
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FIG. 2 illustrates a digital filter forming part 
of the vector estimation system of FIG. 1; 

FIG. 3 illustrates a parameter estimator forming 
5 part of the vector estimation system of FIG. 1; 

FIG. 4 illustrates a smoother module forming part 
of the vector estimation system of FIG. 1; and 

10 FIG. 5 illustrates a speech encoder that includes 

the vector estimation system of FIG. 1. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT OF THE 
INVENTION 

15 

In the drawings, like numerals on different Figs are 
used to indicate like elements throughout. Referring to 
Fig. 1, there is illustrated a vector estimation system 1 
for processing a sequence of input vectors (y 0 to y T ) . The 

20 vector estimation system 1 includes a digital filter 2, a 
parameter estimator 10 and a smoother module 17. The 
digital filter 2 has five inputs and five outputs. The 
five inputs of the digital filter 2 are a filter vector 
input 3, a predictor gain input 4, an observation noise 

25 variance input 5, an OnsetFlag input 2 6 and process noise 
variance input 25. The five outputs of the digital filter 
2 are a current slowly evolving filter estimate output 6, a 
current slowly evolving predicted estimate output 7, a 
previous filter error variance output 9, a previous slowly 

30 evolving filter estimate output 20 and a current predictor 
error variance output 21. 

The parameter estimator 10 has four inputs and three 
outputs. The parameter estimator 10 inputs are an 
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estimator vector input 19 coupled to the vector input 3; a 
previous slowly evolving filter estimate input 13 coupled 
to the previous slowly evolving filter estimate output 20; 
a current predictor error variance input 15 coupled to the 
5 current predictor error variance output 21; and a smoothed 
estimate input 16. The three outputs of the parameter 
estimator 10 are a predictor gain output 11 coupled to the 
predictor gain input 4; an observation noise variance 
output 12 coupled to the observation noise variance input 
10 5; and an OnsetFlag output 22 coupled to the OnsetFlag 
input 26. 

The smoother module 17 has six inputs one being 
coupled to the slowly evolving filter estimate output 6; 

15 one coupled to the slowly evolving predicted estimate 
output 7; one coupled to the previous filter error variance 
output 9; one coupled to the previous slowly evolving 
filter estimate output 20; one coupled to the predictor 
error variance output 21; and one coupled to the predictor 

20 gain output 11. The smoother module 17 also has a smoothed 
estimate output 18 providing an output for the vector 
estimation system 1, the smoothed estimate output 18 is 
coupled to the smoothed estimate input 16 of the parameter 
estimator 10. 

25 

Referring to Fig. 2, the digital filter 2 is a comb 
filter in the form of a Kalman Filter Bank. The digital 
filter 2 comprises a Kalman gain determination module 3 0 
with one input being the observation noise variance input 
30 5. There is also a filtered state estimation module 31 with 
one input being the vector input 3 and another input being 
coupled to an output of the Kalman gain determination 
module 30. An output from the filtered state estimation 
module 31 provides the slowly evolving filter estimate 



16 



output 6 that is coupled to an input of a previous filtered 
state adjustment module 33. Other inputs of the previous 
filtered state adjustment module 33 are provided by the 
OnsetPlag input 2 6 and the vector input 3 via a delay 
5 module 34. 

An output from the previous filtered state adjustment 
module 33 provides the previous slowly evolving filter 
estimate output 2 0 that is coupled to an input of a 

10 predicted state estimation module 35. Another input to the 
predicted state estimation module 3 5 is provided by the 
predictor gain input 4. An output of the predicted state 
estimation module 35 provides the slowly evolving predicted 
estimate output 7 that is coupled to an input of the 

15 filtered state estimation module 31. 

The output from the Kalman gain determination module 
30 is also coupled to an input of a filter variance 
estimation module 32 that has an output coupled to an input 

20 to a previous filter variance adjustment module 36. An 
output from the previous filter variance adjustment module 
3 6 provides the previous filter error variance output 9 
that also provides an input to a predictor variance 
estimation module 37. Other inputs to the predictor 

25 variance estimation module 37 are provided by the predictor 
gain input 4, process noise variance input 25, OnsetFlag 
input 26 and observation noise variance input 5. An output 
from the predictor variance estimation module 37 provides 
the predictor error variance output 21 that is coupled to 

30 inputs of the Kalman gain determination module 30, the 
filter variance estimation module 32 and previous filter 
variance adjustment module 36. Other inputs to the 
previous filter variance adjustment module 3 6 are provided 
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by the predictor gain input 4, the process noise variance 
input 25 and the OnsetFlag input 26. 

As will be apparent to a person skilled in the art, 
5 the characteristics of the digital filter 2 are formalised 
in equations (l)-(6) below. 

At an nth input vector y n (a current input vector) of 
the series of input vectors (y 0 to y T ) received by the 
10 system 1, the previous filtered state adjustment module 33 
provides, at the previous slowly evolving filter estimate 
output 20, a previous filtered estimate value x f n _ x . 

The OnsetFlag input 26 is a binary signal input that 
15 indicates whether or not the beginning of a signal segment 
containing a significant amount of harmonic energy 
(determined by a threshold value) has been detected. If 
OnsetFlag input 2 6 is set to a value that indicates that 
the beginning of such a segment has been detected, then the 
20 previous filtered estimate value x f n _ 1 is set to an a 
previous input vector y nl 

Alternatively, if the OnsetFlag input 6 is set to a 
value that indicates that the beginning of a segment 
25 containing a significant amount of harmonic energy has been 
detected, then the previous filtered estimate value x f n _ 1 
is set to a value of the previous input vector y n-1 that 
was provided at vector input 3 . 

30 For the current input vector y n , the digital filter 2 

provides a current predicted estimate value x p n at the 
predicted estimate output 7 . The current predicted 

estimate value x p n is computed according to: 
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Where oc n is a current predictor gain value for the nth 
input vector y , provided at the predictor gain input 4. 

Once the current predicted estimate value x p n is 
computed, then also for the current input vector y n a 
current filtered estimate value x f ^ n is provided at the 
slowly evolving filter estimate output 6. This current 
filtered estimate value x f n is computed according to: 

X f,n = Vn + K n"^n " X p,n> "~ ~ (2) 

Where K n is a current Kalman gain value for the 
digital filter 2 for the nth current input vector y n . 

The Kalman gain value K n is computed according to: 

K = 2 / (2 + a 2 ) (3) 

n p,n p,n v,n 

Where, Z p n is a current predictor error variance 
value provided at the predictor error variance output 21 

2 

for the nth input vector y n . ; and a v n is a current 
observation noise variance value provided at the 
observation noise variance input 5 also for the nth input 
vector y . 

If the OnsetFlag is set to a value that indicates that 
the beginning of a signal segment containing a significant 
amount of harmonic energy has been detected, then the 
30 current predictor error variance value 2 p n is typically 
set to the observation noise variance value o . This 
results in Equation (3) producing a Kalman gain value K 
equal to 0.5. 
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If OnsetFlag is set to a value that indicates that the 
beginning of a signal segment containing a significant 
amount of harmonic energy has not been detected, then the 
5 current predictor error variance value E p n is computed 
according to: 

S n = a .a .2. n . + a 2 (4) 

p,n n n r,n-J. w 

2 

10 where G w is a process noise variance value provided at 

the process noise variance input 25 and X f n _ 1 is a previous 
filtered error variance value for a previous input vector 

15 If the OnsetFlag is set to a value that indicates that 

the beginning of a signal segment containing a significant 
amount of harmonic energy has not been detected then a 
current filtered error variance value S f n , provided at the 
output of the filter error variance estimation module 32, 

20 is computed according to: 

2 fjn = (1 - K n ) . 2 p|n (5) 



If the OnsetFlag is set to a value that indicates that 
25 the beginning of a signal segment containing a significant 
amount of harmonic energy has been detected, then the 
current filtered error variance value Z f is computed 
according to: 

30 S f n = (Z p ^ n - ct w 2 ) / a n 2 (6) 



Referring to Fig. 3 there is illustrated the parameter 
estimator 10 that typically comprises an initial parameter 
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estimation module 40, an unvoiced speech adjustment module 
41 and a voicing onset adjustment module 42. The initial 
parameter estimation module 40 has four inputs provided by 
the predictor error variance input 15, the slowly evolving 
5 filter estimate input 13, the vector input 19 and the 
smoothed estimate input 16. Outputs of the parameter 
estimation module 40 are coupled to inputs of the unvoiced 
speech adjustment module 41 and a further inputs to module 
41 are provided by the predictor error variance input 15, 

10 the slowly evolving filter estimate input 13 and the vector 
input 19. Outputs of the unvoiced speech adjustment module 
41 are coupled to inputs the voicing onset adjustment 
module 42 and a further input to module 42 is provided by 
the vector input 3 . The voicing onset adjustment module 42 

15 has three outputs providing the predictor gain output 11, 
observation noise variance output 12 and OnsetFlag output 
22. 

The initial parameter estimation module 40 computes 
20 initial estimates of the current predictor gain value a 
and the current observation noise variance value G vn . 
These are determined as follows: 

■ X f,n-l ) 1 (X f,n-l T ' X f f n-1> (7) 

25 

= (y n T . (y n - a^^) /N) - 2p, n (8) 

where N is the number of elements in the current input 
vector y n that is the nth one of the sequence of input 
30 vectors (y 0 to y T ) received at vector input 3; y is the 
transpose of the current input vector y n ; and x f n _ 1 T is the 
transpose of the previous filtered estimate value x f n _ x 
resulting from the previous input vector Y n _ x applied to 
the digital filter system 1. These estimate values are 
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based on the assumption that both the current and at least 
the previous two input vectors represent signal segments 
containing a significant amount of harmonic energy. 

5 It should be noted that for certain applications, the 

current predictor gain value (X n may by a set to a constant 
value of 1. Thus, current observation noise variance value 
o" vn is determined as follows: 

10 O vn 2 = (y n T . (y n - x ftn _ 1 )/N) - Zp, n (8a) 



The unvoiced speech adjustment module 41 determines 

whether the current input vector y n represents a segment of 

15 speech that contains no significant harmonic energy, and if 

so selectively sets the current predictor gain value oc n and 

2 

the current observation noise variance value o V n to 
appropriate values. Preferably, the unvoiced speech 
adjustment unit determines that the current input vector y n 
20 represents a segment of speech that contains no significant 
harmonic energy by detecting whether either of the 
following conditions is true: 

(i) a is less than 0.0, or 

25 (ii) both a n is greater than 1.0 and the initial 

estimation of the observation noise variance 
value is greater than mean squared value of 
elements in the current predicted estimate 
value . 



the 
1.0, 



If either conditions (i) or (ii) hold, then typically 

unvoiced speech adjustment module 41 will set a to 
2 

and re-compute o vn accordingly using Equation (8) . 
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The voicing onset adjustment module 42 determines if 
the current input vector y n represents the second cycle of 
a segment of speech containing a significant amount of 
harmonic energy, and if so adjusts current predictor gain 
5 value a n and the observation noise variance value c V n to 
more appropriate values and sets the OnsetFlag to a value 
indicating that voicing onset has been detected. 

Typically, the voicing onset adjustment module 42 
10 determines that the current input vector y is the second 
cycle of a segment of speech containing a significant 
amount of harmonic energy as follows. An input prediction 
gain, 8, is computed according to: 

15 P^n" • ^n-l 5 / ^n-l -yn-l } (9 > 

An input prediction error variance value, o g , is computed 
according to: 

o- e = y n T . (y n - P-y^) / N (10) 

The voicing onset adjustment unit determines whether both 
20 of the following conditions are true: 

... 2 2 

(m) o" e is less than k^cr^ , wherein k 1 is a 

constant, whose value is typically 0.9. 

(iv) o e divided by the mean squared value of the 

25 elements of the input vector is less than 

k 2 , wherein k 2 is a constant, whose value is 

typically 0.5. 

If both conditions (iii) and (iv) hold, then typically the 

30 voicing onset adjustment unit will set a to 8 and set G 
2 n TO 

to o e . 

Referring to Fig. 4 there is illustrated the smoother 
module 17 that typically comprises series coupled smoothed 
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state estimation modules 5 0 a first stage of which has an 
input receiving the current filtered estimate value x £ n . 
The final stage of the smoothed state estimation module 50 
provides a smoothed estimate value x f at output 18 of a 

5 previous slowly evolving component. The smoother module 17 
also has five sets of series coupled delay modules 
51,52,53,54 and 55 with respective outputs of an jth delay 
module 51,52,53,54 and 55 providing inputs to an j+lth 
smoothed state estimation module 50. 

10 

The smoothed state estimation modules 50 provide 
smoothed estimates X s {n _^ j) for successive values of j 
beginning with j=l. These estimates are computed according 
to: 

15 

X s,(n-j, = -f^n-j) + C -<*s,(n- j + l) " X p,(n- j+ l>> —-(11) 

wherein 

c = ^f.-rWWiKi,) —-(12) 

20 and 

X s,n = X f,n (13) 

From the above it will be apparent that the purpose of 
25 the smoother module 17 is to provide an estimate X s (n _ 3) of 
the slowly evolving component of an input vector y n _. based 
upon input vectors up to and including y n . The smoother 
module 17 thus uses current data to estimate a past slowly 
evolving component value, in contrast to the digital filter 
30 2, which uses current data to estimate a current slowly 
evolving component value. 
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In use, the vector estimation system 1 receives the 
sequence of input vectors y 0 to y T that are each comprising 
N elements. Each of the input vectors y Q to y T contains a 
sampled period of a presumed quasiperiodic signal. This 
5 sampled signal is typically time warped to allow for 
variations of quasiperiodic periods, so that each input 
vector contains the same number of elements, as will be 
apparent to a person skilled in the art. Alternatively, 
consecutive input vectors y Q to y T may have elements added 

10 to them or removed from them, again so that the resulting 
number of elements in each is the same. For an nth 
iteration, a input vector y n will be applied to vector 
input 3 and estimator vector input 19. The digital filter 
2 processes this input vector y n resulting in the slowly 

15 evolving filter estimate output 6 providing, to input 13, 
the previous filtered estimate value x f n _ 1 of a slowly 
evolving component of sequence of vectors y Q to y T . 

The parameter estimator 10 processes the previous 
20 filtered estimate value x f n _ 1 and current input vector y n 
to provide a current current predictor gain value oc n at 
predictor output 11. The current predictor gain value an 
is thereby applied to input 4 of the digital filter 2 for 
controlling the gain thereof during filtering of input 
25 vector y . The parameter estimator 10 determines the 
current predictor gain value' OL n by the calculation stated 
in equation (7) . 

As will be apparent to a person skilled in the art, at 
30 initialisation (i.e. the first sample time when n is 0 
therefore input vector y 0 is applied to digital filter 
system 1), there will be no previous filtered estimate 
value x f n _ 1 . Accordingly, although there are many ways to 
allocate a value for this previous filtered estimate value 
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x f n _ 1 , the present invention preferably assigns the 
previous filtered estimate value x f n _ 1 with the same 
element values as input vector y Q . 

5 Referring to Fig. 5, the vector estimation system 1 

can advantageously included in a speech encoder 60. The 
speech, encoder 60 includes a signal normalization module 61 
with an input for receiving a speech signal. A signal 
vector output 62 of the signal normalization module 61 is 

10 coupled to filter vector input 3 and another output is 
coupled to the process noise variance input 25 of the 
vector estimation system 1. The signal vector output 62 is 
also coupled to an input of an adder module 63 and another 
input of the adder module 63 is coupled to the smoothed 

15 estimate output 18 of the vector estimation system 1. An 
output from the adder module 63 is coupled to an input of a 
rapidly evolving component encoder 64 and there is also a 
slowly evolving component encoder 65 having an input 
coupled to the smoothed estimate output 18. The speech 

20 encoder 60 has three outputs from the signal normalization 
module 61, for coupling to a speech decoder, these outputs 
being a spectral envelope output 66, an energy parameters 
output 67 and a period parameters output 68. The speech 
encoder 60 also has a slowly evolving component output 69 

25 from the evolving component encoder 65 and a rapidly 
evolving component output 70 from the rapidly evolving 
component encoder 64 . 

In operation, the speech encoder 60 firstly normalizes 
30 a speech signal with respect to its spectral envelope, 
energy and period. The normalisation process involves 
estimating parameters that describe the spectral envelope, 
energy and period of the input signal and these parameters 
are typically transmitted to a speech decoder at outputs 
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66,67,68. Process noise variance provide at the process 
noise variance input 2 5 is typically used to control the 
vector estimation system 1. The normalisation process 
produces the sequence of input vectors (y 0 to y T ) for the 
5 vector estimation system 1. The sequence of input vectors 
(y 0 to y T ) are a sequence of fixed length vectors 
representing sampled consecutive cycles of the normalised 
waveform. These vectors (y 0 to y T ) are applied to the 
filter vector input 3 of the vector estimation system 1, 

10 which generates a slowly evolving component at the smoothed 
estimate output 18. By subtracting this slowly evolving 
component from the input vectors (y 0 to y T ) , a rapidly 
evolving, or noise-like component is produced and provided 
to the rapidly evolving component encoder 64. The slowly 

15 evolving and rapidly evolving components are encoded 
respectively by the slowly and rapidly evolving component 
encoders 64, 65. The encoders 64,65 use appropriate 
methods known in the art to produce parameters at 
respective outputs 70,69 which are transmitted to a speech 

20 decoder. 

Advantageously, as can be seen with reference to 
equation (7) , the present invention provides for the vector 
estimation system 1 to receive the current input vector y n 

25 that is one of a the sequence of input vectors y 0 to y T . The 
parameter estimator 10 provides the current predictor gain 
value a n at the predictor gain output 11 thereby modifying 
the current filtered estimate value x f n at the slowly 
evolving filter estimate output 6. In this regard, the 

30 current predictor gain value a is dependent upon the 
previous filtered estimate value x f n _ 1 and the current 
input vector y n . As will be apparent to a person skilled 
in the art, the parameter estimator 10 determines the 
current predictor gain value CC n from both a sequence of 
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previous input vectors Y u _ 1 to Y 0 an -d a sequence of 
corresponding filtered estimate values provided by 
successive preceding values of the previous filtered 
estimate value x f 0 to x f n _ 1 • 

5 

With reference to equation (8) , the present invention 
also advantageously provides for the vector estimation 
system 1 to receive the current input vector y n The 
parameter estimator 10 provides the current observation 

10 noise variance value <7 v n 2 at the observation noise variance 
output 12 thereby modifying current filtered estimate value 
x f n at the slowly evolving filter estimate output 6. The 
current observation noise variance value 0" v n 2 is dependent 
upon the current input vector y n , the current predictor 

15 gain value a n (that can be set to a constant value of 1 for 
certain applications) , the current predictor error variance 
value £p,n/ and said previous filtered estimate value x f n _ 1 . 

The detailed description provides a preferred 
20 exemplary embodiment only, and is not intended to limit the 
scope, applicability, or configuration of the invention. 
Rather, the detailed description of the preferred exemplary 
embodiment provides those skilled in the art with an 
enabling description for implementing a preferred exemplary 
25 embodiment of the invention. It should be understood that 
various changes may be made in the function and arrangement 
of elements without departing from the spirit and scope of 
the invention as set forth in the appended claims. 



